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This is a set of four lectures presented at the Theoretical Advanced Study 
Institute (TASI-09) in June 2009. I provide an introduction to experiments 

designed to detect WIMP dark matter directly, focussing on building intuitive 
understanding of the characteristics of potential WIMP signals and the ex- 
perimental techniques. After deriving the characteristics of potential signals 
Q in direct-detection experiments for standard WIMP models, I summarize the 

^ general experimental methods shared by most direct-detection experiments 

and review the advantages, challenges, and status of such searches (as of late 
2009). Experiments axe already probing SUSY models, with best limits on the 
spin-independent coupling below 10~'^pb. 

Keywords: WIMPs, cold dark matter, SUSY, dark matter halo 

> 

in 

1. Introduction 

<N 

lO A variety of astrophysical observations-^'^ indicate that 83% of the matter in 

^-H the Universe is nonbaryonic and dark, presumably in the form of elementary 

particles produced in the early Universe. Because no such particles have yet 
been identified in particle accelerators, these observations require new fun- 
• • damental particle physics. Weakly Interacting Massive Particles (WIMPs^) 

. ^ form a particularly interesting generic class of new-particle candidates be- 

cause they naturally provide about the inferred amount of this nonbaryonic 
^ dark matter,^ a result dubbed the "WIMP miracle." 

WIMPs would be produced thermally in the early Universe. Because 
they interact only weakly, their annihilation rate would become insignifi- 
cant as the Universe expands, thus "freezing out" a relic abundance of the 
particles (for a pedagogical discussion, see e.g. Refs. 5-8). The expected 
WIMP density would be the same as that of the nonbaryonic dark matter 
if the WIMP velocity- averaged annihilation cross section is ^Ipb, so that 
the WIMP mass is -lOOGeV/c^. 
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Remarkably, extensions to the Standard Model motivated entirely by 
particle physics predict particles with the same cross section and mass as 
these dark-matter WIMPs. Detection of the W and Z bosons indicates 
that the electroweak symmetry is spontaneously broken at a scale of ~ 
100 GcV/c^. Whatever physics solves the hierarchy problem associated with 
this symmetry breaking - be it supersymmetry,^'®^^^ extra dimensions/^""'^^ 
or something else - gives rise to additional particles. If an appropriate 
(often independently motivated) discrete symmetry exists, the lightest such 
particle is stable. This particle is then weakly interacting, massive, and 
stable - it is a WIMP. Thus, particle theorists are "almost justified in 
saying that the problem of electroweak symmetry breaking predicts the 
existence of WIMP dark matter". ^"^ 

Although the argument for WIMP dark matter is generic, supersymme- 
try dominates the discussion as a particularly well-motivated model. Su- 
pcrsymmctry interactions arise in theories of quantum gravity, stabilize the 
Higgs mass hierarchy problem, predict the observed value of sin^ and 
over a broad range of parameter space predict cosmologically significant 
relic WIMP densities. 

WIMPs can potentially be detected by three complementary methods. 
They may be produced and detected (indirectly) at accelerators such as 
the Large Hadron Collider (see e.g. Ref. 18). Relic WIMPs may be de- 
tected indirectly when they clump in massive astrophysical objects, increas- 
ing their annihilation rate enough that their annihilation products may be 
detectable. Many potential (or suggested^^"^^) indirect signals are 
ambiguous, with alternate astrophysical explanations (see Refs. 7,29 and 
references therein). Some potential indirect signals, however, would be com- 
pelling. Annihilation in the Sun or Earth would produce higher-energy neu- 
trinos than any other known process. These neutrinos could be observed in 
neutrino telescopes such as IceCube^^'^^ or ANTARES.-'^^.s.s Y^ithei FERMI 
or ground-based air Cerenkov telescopes may detect distinctive gamma-ray 
features from the galactic center or from sub-halos.'^'*"^^ 

Relic WIMPs may also be detected directly when they scatter off nu- 
clei in terrestrial detectors. '^^^^^ This article offers an introduction to these 
direct-detection experiments. Section 2 includes derivations and explana- 
tions of the characteristics of potential signals in direct-detection experi- 
ments for standard WIMP models. Section 3 summarizes the general exper- 
imental methods shared by most direct-detection experiments and discusses 
particular experiments briefly, emphasizing the relative advantages and dif- 
ferent challenges and capabilities of the various approaches. 
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2. WIMP-nucleus elcistic scattering: from model to signal 

Understanding experiments designed for direct detection of dark matter 
begins with the observables of potential signals. In this section we consider 
the observables of any model that predicts standard WIMP-nucleus elastic 
scattering (see Neil Weiner's contribution to these proceedings^ for discus- 
sion of more speculative models with non-standard scattering). Following 
the reviews by Lewin and Smith, and Jungman, Kamionkowski and Gri- 
est,^ this section derives how the observed WIMP interaction rate depends 
on energy, target, time, and direction. 

2.1. Spin-independent and spin- dependent cross sections 

Using Fermi's Golden Rule, we can divide the energy dependence of the dif- 
ferential WIMP-nucleon cross section into a term aowN that is independent 

of the momentum transfer and a term F'^(q) (known as the form factor) 
containing the entire dependence on the momentum transfer q: 

dawti{q) _ 1 I . .,2 _ o"owNf^(g) 

dq^ ~ TTV^^^ ~ An\v^ • ^ ' 

Here, v is the velocity of the WIMP in the lab frame, and the WIMP- 
nucleus reduced mass = M^MA/{My^ + Ma) in terms of the WIMP 
mass and the mass Ma of a target nucleus of atomic mass A. Since the 
WIMPs are nonrelativistic, the zero-momentum cross section for a WIMP 
of arbitrary spin and general Lorentz- invariant WIMP-nucleus cross section 
may be written in terms of a spin-independent (mostly scalar) and a spin- 
dependent (mostly axial vector) term: 

aowN = ^ [ZU + Z)Uf + ^ (ap(5p) + a^{S^)f . 

TT n J 

(2) 

The proof of this claim makes a good exercise for the reader; solution may be 
found in Ref. 42. Here /p and /„ (op and an) are effective spin-independent 

(spin-dependent) couplings of the WIMP to the proton and neutron, respec- 
tively. Together with the WIMP mass, M-^, these parameters contain all the 
particle physics information of the model under consideration. The other 
parameters describe the target material: its atomic number Z, total nuclear 
spin J, and the expectation values of the proton and neutron spins within 
the nucleus {Sp^n) = {^\Sp,n\N)- For free nucleons, {Sp) = (S'„)= 0.5. Ta- 
ble 1 from Ref. 43 lists values of {Sp) and {Sn) for materials commonly 
used for dark matter searches, although some are subject to significant 
nuclear-physics uncertainties. 
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Table 1. Values of the atomic number Z, the total nuclear spin J, and the expec- 
tation values of the proton and neutron spins within the nucleus {5p,n} for various 
nuclei with odd numbers of protons or neutrons, leading to the relative sensitivities 
to spin-dependent interactions shown, from Refs. 5,43 and the references contained 
therein. 
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For many models, /p ~ /n, so the spin- independent WIMP-nucleus cross 
section 

TT 

The dependence of this cross section on the target material may be factored 
out by rewriting this result as 

"■qwn.si = o'si^^^, (4) 

where /in is the reduced mass of the WIMP-nucleon system, and the (target- 
independent) spin-independent cross section of a WIMP on a single nucleon 

crsi = . (5) 

TT 

This WIMP-nucleon cross section asi may be used to compare experimen- 
tal results to theory and to each other. A given model predicts particular 
combinations of usi and My/, different experiments produce limits on crgi as 
functions of by translating limits on the WIMP-nucleus cross-section 
to limits on crsi using equation 4. The dependence on /i^A^ in eqn. 4 indi- 
cates the advantage of experiments using relatively heavy target materials 
(but see the effects of the form factor in Sec. 2.3). For a SOGcV/c^ WIMP 
incident on a target with A = 50, = 625, so the spin-independent 
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WIMP-nucleus cross section is larger than the WIMP-nucleon cross section 
by a factor > 10''. 

The situation for spin-dependent interactions is quite different.*^ First 
of all, contributions from the spin-dependent proton and neutron couplings 
often cancel, so it is important to quote limits on the spin-dependent in- 
teraction on neutrons separately from that on protons, each under the as- 
sumption that the other interaction is negligible. Furthermore, while the 
coherent interaction on the nucleus results in a spin-independent cross sec- 
tion that scales with since the contribution of each nucleon adds inside 
the matrix element, the spin-dependent contributions of nucleons with op- 
posite spins cancel, so that the total spin-dependent cross section depends 
on the net spin of the nucleus. As shown in Table 1, nuclei with even num- 
bers of protons have nearly no net proton spin and essentially no sensitivity 
to spin-dependent interactions on protons, and nuclei with even numbers 
of neutrons similarly have almost no sensitivity to spin-dependent interac- 
tions on neutrons. Argon, with even numbers of protons and neutrons for 
all significant isotopes, is thus insensitive to spin-dependent interactions. 
Many materials used as WIMP targets (Ge, Si, Xe) have even numbers of 
protons and hence are insensitive to spin-dependent interactions on pro- 
tons; only some isotopes of these targets (and hence only a fraction of the 
detector's active mass) have sensitivity to spin-dependent interactions on 
neutrons. Typically, sensitivity to spin-dependent interactions on protons 
requires alternate target materials, often resulting in worse backgrounds or 
background rejection and lower sensitivity to spin-independent interactions. 
The relative sensitivity of a material to spin-dependent interactions is sum- 
marized by its "scaling factors" 4(S'p)2(J + 1)/3J and 4(S'„)2(J + 1)/3J, 
which are listed in Table 1. As with spin-independent limits, experimenters 
quote limits on target-independent quantities: the spin-dependent WIMP- 
proton cross section asup = 24G'|/ipap/7r and the spin-dependent WIMP- 
neutron cross section ctsdii = 24G|/iJ;a^/7r. 

The lack of benefit from the coherent interaction for spin-dependent 
interactions results in most models being more accessible experimentally 
via their spin-independent interactions than by their spin-dependent inter- 
actions. As shown in Fig. 1 and Fig. 2, current experiments are already 
constraining MSSM models based on their spin-independent couplings, but 
none is yet sensitive enough to constrain such models based on their spin- 
dependent couplings, despite the fact that spin-dependent couplings are 
typically ~ 3000 x larger than spin-independent couplings. 
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Fig. 1. Upper limits on the spin-independent WIMP-nuclcon coupling crgi under 
the standard assumptions about the Galactic halo described in the text. Most sensi- 
tive limits are from cryogenic experiments (soUd) CDMS** (black), EDELWEISS-II'''5 
(medium gray), and CRESST*^ (light gray), and two-phase noble experiments (dashed) 
XENONIO*'' (black), ZEPLIN-III^s (medium gray), and WArP^s (light gray). Current 
experiments already exclude part of the parameter space of MSSM models (shaded).®" 
Figure made using the Dark Matter Limit Plotter. 

2.2. The WIMP recoil energy spectrum 

It is illuminating to calculate the energy spectrum for the case of zero 
momentum-transfer (i.e. taking = 1). Furthermore, simply multiplying 
this spectrum by the energy dependence of F^{q), rather than including 
the form factor F"^ within the kinematic integral to follow, is convenient 
and usually adequate. 

The energy spectrum arises due to the familiar kinematics of elastic 
scattering. In the center-of-momentum frame, the WIMP scatters off a nu- 
cleus through an angle 6, with cos^ uniformly distributed between —1 and 
1 for the isotropic scattering that occurs with zero-momentum transfer. 
If the WIMP's initial energy in the lab frame Ei = the nucleus 

recoils with energy 



_ (1-cos^) 

— i^iT ^ 



(6) 
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Fig. 2. Upper limits on the spin-dependent WIMP-neutron coupling o"sd,j (left) and 
the spin-dependent WIMP-proton coupling asiip (right) under the standard assump- 
tions about the Galactic halo described in the text. The most sensitive limits on trsDn 
are from the same experiments shown in Fig. 1 (with the same linetypcs): XENONIO''"'^ 
(black dashes), ZEPLIN-III53 (medium gray dashes), and CDMS*'' (black solid). Note 
ZEPLIN-III limits were calculated with a scaling factor 2x smaller than that used for 
XENONIO. Due to the low intrinsic sensitivity of leading (Xe and Ge) experiments 
to spin-dependent interactions on protons, the most sensitive limits on (TSDp are from 
experiments with only modest sensitivity to spin-independent interactions: PICASSO'"'^ 
(6-sided stars), COUPP^"' (5-pointed stars), KIMS"'^ (circles), and NAIAD'" (x). Limits 
from indirect search experiments SuperKamiokande''** (points) and IceCube'^-' (dotted) 
make additional assumptions about branching fractions to neutrinos. Current exper- 
iments do not exclude any part of the parameter space of the same MSSM models 
(shaded)^*^ shown in Fig. 1, despite the fact that the predicted spin-dependent cross sec- 
tions axe ~ 3000 X larger than the spin-independent ones. Figure made using the Dark 
Matter Limit Plotter.^^ 



(in the lab frame), where 

M^Ma (Mx + MAf 

is a dimcnsionlcss parameter related to the reduced mass fiA- Note that 
r < 1, with r = 1 only if M-^ = Ma- For this isotropic scattering, the 
recoil energy is therefore uniformly distributed between 0-Eir. As shown in 
Fig. 3, the differential contribution to the differential rate for a given initial 
WIMP energy 



(8) 
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Fig. 3. Plot showing schematically how contributions to the diflferential rate dR/dE^ 
for different values of the initial WIMP energy add. We define total differential rate 
dR/dE^ of WIMPs with initial energy to be dR{Ei). For a WIMP initial energy 
E{, the recoil energy E^i is uniformly distributed between 0— -Eir, so dR(Ei) contributes 

equally to the rates of all recoils between E^r, as depicted by the shaded area in the 
figure. The contribution to the differential rate at a given recoil energy (the height of 
the shaded area in the figure) is simply the area divided by the length, or dR(E{) / (E^r) . 
The total differential rate can then be found by summing all the boxes, i.e. integrating 
dR{Ei)/{Eir) for all Ei. 



The maximum initial WIMP energy may be taken as infinity as an initial ap- 
proximation, or more accurately may be based upon the Galactic escape ve- 
locity, Vesc- -Emax = -^x"csc/2- To cause a recoil of energy Er, the minimum 
initial WIMP energy E^min = Eji/r (for head-on scattering, with = tt), 
and the minimum WIMP velocity Vmin = ^/2E^^~JM^ = ^/2E^J{rM^. 

To determine the rate of WIMP-nucleus scattering, it is helpful to imag- 
ine the motion of the target nucleus relative to WIMPs with velocity v in 
the lab frame. In time dt, each nucleus interacts with any WIMP inside 
a volume dV = avdt, where a is the WIMP-nucleus cross section. The 
number of WIMPs inside the volume moving with velocity v 



where the local WIMP number density uq = p^/M^, where is the mass 
density of WIMPs in the galaxy, estimated from studies of Galactic dynam- 



so 





dN = nof{v + v-E)(jvdt, 



(10) 
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ics to be about 0.3 GcV/ (c^ cm"^) (with wide systematic uncertainties^^"^^). 
Note that the number density based on = 0.3 GeV/ (c^ cm^) is really an 
upper limit, since Galactic dark matter may include species other than 
WIMPs. We use the fact that the velocity Vg of the WIMP in the galaxy 
is the vector sum of the WIMP velocity with respect to the Earth v and 
the velocity ve, of the Earth with respect to the Galaxy. We assume that 
the WIMPs' velocities in the frame of the Galaxy follow the Max;wellian 
distribution: 

fiv + VE) = ^ (11) 

where Vq = 220 ± 20km/s is the local circular velocity, ^'^ and A; is a nor- 
malization factor. This simple distribution is not expected to be especially 
accurate, but it provides a useful standard. Sec e.g. Refs. 64,65 for discus- 
sion of alternatives that are likely more accurate, and see Section 2.4 below 
for discussion of the impact of uncertainties on a WIMP discovery. 

The differential interaction rate per kilogram of detector is then the 
product of the number of interactions per nucleon with the number of nuclei 
per kg of material: 

dR = ^nof{v + VE)av(fv (12) 

where A^o is Avogadro's number, so that Nq/A is the number of nuclei per 
kilogram of material. 

It is instructive (and reasonably accurate, as shown below) to consider 
the simplified case ignoring the Earth velocity and the Galaxy's escape 
velocity {i.e. setting ve = 0, Vesc = oo), for which the integral is trivial. 
After setting 

Ro = -^—r-noav, (13) 

VTT A 

we get 

~ Eor ' ^ ' 

where Eq = is the most probable WIMP incident energy. The 

mean recoil energy is easily seen: (-Er) = Egr. Since r < 1, the mean recoil 
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{En) = Eor = -M^vl « 15 keV. (17) 



10 

energy {E-r) = Eq only if the WIMP mass is equal to the mass of the target 
nucleus; (-Er) < Eq both for smaller and for larger WIMP masses. As an 
example, since vq « 220km/s w (0.75 x 10-^)c, = Ma = 50GeV/c^ 
would result in 

1 

A different target mass would result in even lower {En)- This low energy 

sets the; first challciiige for direct detection experiments they must have 
low energy thresholds, much lower than past solar neutrino experiments for 
example. 

From the exponential form of the approximate energy spectrum, we see 
that Rq is the total WIMP rate. If we plug known numerical values into 

equation 13, we find 

D ^ (^OWN Px . 1 -Ij -1 /10^ 

■Ro ~ , ^ / ox ^ — , events kg day . (18) 

Mx(GeV/c2) Ipb 0.4GeV/cm^ ^ ' 

A 50GeV/c^ WIMP with a WIMP-nucleus cross section (Jqwn = 1 pb (so 
that the spin-independent WIMP-micleon cross section ctsi ^ 10^^ pb, or 
the spin-dependent WIMP-nucleon cross section cisDp.n ~ 10""^ pb) results 
in about 10 events/ (kg day). Since the energy spectrum is a falling expo- 
nential, a low energy threshold is critical to detect most of these events; the 
fraction of events above an energy threshold i?th is e"^"'/^"'". 

The dependence of the energy spectrum on the WIMP mass may be 
seen easily from equation 16. The mean recoil energy 

(En) — rEn oc ^ oc < .n,^ ,^ • (19) 

^ ' + \ constant if > Ma ^ ' 

Heavy WIMPs all yield about the same energy spectrum. This result holds 
for calculations made including the Earth velocity. Galaxy escape velocity, 
and nuclear form factor, as shown in Fig. 4. 

WIMPs with velocities above the Galaxy's escape velocity are likely to 
have already escaped. The finite escape velocity (~ 540km/s « 2 x 10~"^c 
according to the RAVE survey^^) alters the recoil spectrum slightly and 
produces a cut-off at 

En.e^ = ^rMy^^^ « 100 keV. (20) 

The cutoff energy has the same dependence on the WIMP mass as the mean 
recoil energy (see equation 19) since 

E^^='^{En)'=^6{En). (21) 
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Fig. 4. Expected interaction rate on Ge as a function of recoil energy for different 
WIMP masses M^j, (as shown in the legend in units of GeV/c^), for most probable 
WIMP velocity tig = 270km/s, on the high end of possible values. The spectra of lower- 
mass WIMPs are softer and cut off at lower energy due to the Galactic escape velocity. 
The spectra of WIMPs heavier than the target nucleus are nearly identical. The rate is 
lOx larger for a lOOGeV/c^ WIMP than for a lOOOGeV/c^ WIMP since there would be 
lOx more lOOGeV/c^ WIMPs than lOOOGeV/c^ WIMPs. Higher-mass spectra deviate 
from straight lines due to non-unity form factor (see Section 2.3). 



Hence, higher-mass WIMPs produce recoils that are easier to detect and 
have cutoff energies so high as to usually be negligible. The cutoff energy, 
however, is significant for low-mass WIMPs (see Fig. 4); experiments will 
have no sensitivity at all to WIMPs of low enough masses due to the cutoff. 
Since the Galactic escape velocity is not known especially well, caution 
should be taken when drawing conclusions that may be sensitive to the the 
number of WIMPs with velocities at or near the assumed cutoff energy. It 
should also be noted that for historical reasons (in order to quote limits 
using the same assumptions as previous experiments) most experimenters 
routinely assume the "standard halo" described in Ref. 67, which uses a 
value Vcsc = 650km/s, somewhat above the 90% upper limit quoted in the 
more recent RAVE survey. 

The dependence of the energy spectrum on target mass Mj^ (ignoring 
the form factor F'^) is entirely through the r parameter in Eqn. 19 or 
Eqn. 20. For a given WIMP mass, the cut-off energy i?niax and the mean 
recoil energy are largest for targets whose masses most closely match the 
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Fig. 5. Left: Comparison of the spectrum using the exponential approximation (eqn. 23 
with no cut off energy -Emax) to the fuU calculation (eqn. 22) for a WIMP with mass 
M = 38GeV (solid overlapping curves). The residuals are shown (dashed) after scal- 
ing by 100. At higher energies, the fractional deviations become larger but the absolute 
deviations remain small. The shape of the residuals > SOkeV is from applying the cut- 
off energy Emux in the full calculation but not to the exponential approximation. The 
approximation is slightly less accurate for lower-mass WIMPs. Right: Dependence of 
spectral fit parameter C2 on WIMP mass for a Ge detector with a 10-keV threshold. 
Use of the advised approximate value, C2 = 0.561 (dashed), produces not more than a 
50% error for all WIMP masses. 



WIMP mass. As shown below, including the form factor makes the energy 
spectra of more massive targets softer. 

The full calculation of the energy spectrum for WIMP-nucleus elas- 
tic scattering including the effects of both escape velocity and the earth's 
velocity is left as an exercise for the reader (for an almost complete so- 
lution see Lewin and Smith**^). The result for recoil energies such that 

(^^r) < Uosc - VE is 
dR Rq ( Vq^/tT r / Winin + WE 



R 



Eor I Av 



rf 



— erf 



I'E 



c/"Q 



while for Vcsc 
become^®'^^ 



(22) 

ve < v„iin(E^) < Vesc + ve the terms inside curly braces 



4fE 



erf 



Wcsc 



erf 



WE 



Wo 



WE ■ 



2WE 



This energy spectrum is reasonably approximated by another falling expo- 
nential: 



dR 

dE^. 



[Eb 



Eor 



(23) 
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as shown in Figure 5. Here ci « 0.75 and C2 ~ 0.56, although both depend 
on WIMP and target masses, day of year, and the energy range of interest."'^ 
For example, the right panel of Fig. 5 shows the dependence of C2 on the 
WIMP mass for a Ge target. 

Since Ci/c2 = 1.3, the Earth's motion increases the interaction rate by 
~30% in addition to making the energy spectrum harder, as should be 
expected from analogy to a car's driving through the rain resulting both in 
more raindrops hitting the front windshield and in the drops hitting with 
more force on average. Despite the wide use of this analogy, it must be 
noted that the effect of Earth moving through the WIMP rain is not nearly 
as pronounced as when one drives a car through rain, since the WIMPs are 
moving with velocities comparable to the Earth's. 

The Earth's velocity in the Galaxy of course is not constant over a year, 
but varies due to the small velocity of the Earth around the Sun. As a 
function of the day of the year t, 



with maximum occurring at t — 152.5 days, or June 2. From equation 22, 
one can show that dR/dvE ~ R/2ve, so the 6% annual modulation in the 
Earth velocity from equation 24 causes about a 3% annual modulation in 
the total WIMP interaction rate. Note that this result is true only when 
considering all interactions, even those down to zero recoil energy. The 
interaction rate above a (non-zero) experimental threshold energy can be 
as big as 7%. 

The motion of the Earth in the Galaxy towards the constellation Cygmis 
makes the WIMP flux in the lab frame sharply peaked, resulting in a higher 
rate of recoils from the direction of Cygnus. For the standard halo model 
considered here (and neglecting the escape velocity),™ 



where ^ is the recoil angle in the laboratory relative to the direction of 
Cygnus. Since the Earth speed is comparable to the mean WIMP speed, 
the rate in the forward direction is roughly an order of magnitude larger 
than the rate in the backward direction.™ 




(24) 




(25) 
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Fig. 6. Assumed density of scattering centers for spin-independent interactions, as pro- 
posed by Helm.^^ Density is constant within the nuclear radius rn then decreases to 
zero over a skin thickness s (the related 10%— 90% thickness t is shown in this diagram). 
The Fourier transform of this distribution yields the Woods-Saxon form factor used for 
spin-independent scattering. 



2.3. Nuclear Form Factors 

Under the approximation of plane-wave (Born) scattering, 

M{q) = fnA [ d^xp{x)e''l-\ (26) 



We may identify the momentum-dependent part of this interaction, the 
form factor 

F{q) = j d''xp{x)e^?-\ (27) 

as the Fourier transform of the scattering site positions. For spin- 
independent interactions, a good approximation^^ is the Woods-Saxon form 
factor 

pl^^-^ ^ 3 [sin(grn) - qr^ cos(grn)] ^-(^s?/2^ ^28) 

which is the Fourier transform of a sohd sphere of radius with a skin 
thickness s, as shown in Figure 6. In practice, Lewin and Smith'*^ recom- 
mend values of s = 0.9 fm and 

rl ^ (1.23^1/3 _ Q_60 fm^^ + ^(0.527r fm)^ - 5s^ (29) 



January 28, 2011 1:7 WSPC - Proceedings Trim Size: 9in x 6in SchneeTASIPDF 



15 




Recoil Energy (keV) 

Fig. 7. Spin-independent form factors as a function of recoil energy for targets of 6 
atomic masses A. Prom top to bottom on plot, materials are Ne {A = 20; F or Na are 
similar), Si (A = 28), Ar {A = 40), Ge (A = 73), Xe (dashed, A = 131, I is similar), and 
W (A = 183). 



For spin-dependent interactions, the situation is more complicated. A 
first approximation starts with a thin sheU of valence nucleons, 

= (30) 

n 

but must be extended with detailed nuclear-physics calculations.^^ 

In cither case, F(q) < 1 when the de Broglie wavelength X < Vn and 
the WIMP ceases to interact coherently with the entire nucleus. Since the 
nuclear radius ~ A^/^ fm, this criterion may be rewritten 

, h he 197MeVfm 

q V2Mj^E^ ^2A£R(keV) 

Hence, coherence is lost when 

2 X 10^ 

> keV ~ 100 keV. (32) 

The strong dependence on A indicates that coherence is lost much earlier 
for high-^ targets, as shown in Fig. 7. This loss of coherence significantly 
reduces the advantage of using particularly heavy target materials; practi- 
cally speaking use of materials heavier than Ge yield only modest increases 
in overall rate, far short of the increase that would occur without the 
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20 40 60 80 20 40 60 

Recoil Energy (keV) Threshold Energy (keV) 



Fig. 8. Spin-independent interaction rates (per detector exposure) as a function of 
recoil energy for a WIMP on targets of 6 atomic masses A. From top to bottom on each 
plot, materials are W (A = 183), Xe (dashed, A = 131, I is similar), Ge {A = 73), 
Ar (A = 40), Si {A = 28), and Ne {A = 20, F or Na are similar). Left: Differential 
rate for a 60-GeV/c^ WIMP. High-A materials have a higher rate at low energies, since 
the rate oc fi^A'^, but loss of coherence greatly decreases the rate in these materials 
at high energies. As A increases towards M^, the mean energy and cutoff energy both 
increase due to kinematics, while loss of coherence offsets the increase in the mean 
energy. As A increases past M^, the energy spectrum becomes softer and the cutoff 
energy decreases. Right: Integral rate above the energy threshold indicated for a 100- 
GeV/c^ WIMP. Although energy thresholds vary from experiment to experiment, typical 
energy thresholds for each material are indicated by + signs on each curve. With these 
thresholds, the 100-GeV/c^ WIMP would produce the highest signal rate in Xe, with 
rates in W and Ge about 40% lower. I follows about the same curve as Xe, typically with 
a 3x higher threshold and half the rate. Rates in Si are ~ 9x lower than in Xe, rates in 
Ar are ~ 14x lower, and rates in Ne (or Na or F with this threshold) are ~ 100 x lower. 

loss of coherence. Since the loss of coherence makes these high-^ targets 
intrinsically insensitive to high-energy depositions, it is particularly criti- 
cal that experiments with high-^ materials achieve low energy thresholds. 
Figure 8 shows the relative rates for the same WIMP in several different 
targets. 

2.4. Implications of a detection 

Because the spin-independent, proton-spin-dependent, and neutron-spin- 
dependent form factors are different for a given target, it is possible in 

principle to distinguish the type of interaction by the energy spectrum on a 
single target isotope. Differences are insignificant for low-mass WIMPs since 
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Fig. 9. Expected energy spectra for a 20-GeV/c^ (solid), 60-GeV/c^ (dashed), or 500- 
GeV/c^ (dotted) WIMP interacting on Ge via the neutron-spin-dependent (lowest), spin- 
independent (middle), or proton-spin-dependent (top) interaction. For relatively massive 
WIMPs for which the loss of coherence is significant, a high-statistics detection could 
identify the type of interaction via the spectral shape. Figure provided by J. Filippini. 



all form factors are essentially unity. However, as shown in Fig. 9, differences 
for high-mass WIMPs may be significant with a sufficiently high-statistics 
detection (particularly if the WIMP mass is known independently). 

Detections with several target materials would reveal the type of inter- 
action more clearly. Comparing the interaction rate in different materials 
would indicate the interaction type since the rates scale differently in each 
(Figure 8 shows the scaling for spin-independent interactions, while Table 1 
shows the material-dependent scaling factors for spin-dependent interac- 
tions). Furthermore, detection with different target materials provides a 
useful confirmation of the detection, especially if a consistent WIMP mass 
is determined from each. 

Measurement of the WIMP recoil spectrum would provide constraints 
on the WIMP mass, as can be seen from Fig. 4. However, since heavy 
WIMPs all yield about the same energy spectrum, as shown by Eqn. 19, 
detection of a heavy WIMP would provide only weak constraints on its 
mass, other than that it must be relatively heavy. Table 2 lists how well 
a WIMP mass may be determined by a detection for a spin-independent 
interaction on Ge if the WIMP velocity distribution is known. The uncer- 
tainties on the type of interaction and WIMP velocity distribution would 
contribute additional uncertainty on an inferred WIMP mass. As shown 
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Table 2. Projected limits (at the 99% confidence level) 
on WIMP mass for a GO-GeV/c^ WIMP on Ge based on 

statistical uncertainties only.''^ Also listed is the minimum 
mass Mjnin for which there is no upper limit on the WIMP 
mass. 



Events 


Lower 


Upper 




Detected 


Limit 


Limit 




10 


30 GeV/c2 


none 


50 GeV/c2 


100 


45 GeV/c2 


101 GeV/c^ 


100 GeV/c2 


1000 


55 GeV/c2 


69 GeV/c2 


250 GeV/c2 



in Fig. 9, the different form factors may help a lower-mass WIMP with 
one predominant interaction produce a spectrum more similar to that of a 
higher-mass WIMP. For a Maxwellian distribution, the uncertainty on Vq 
translates into uncertainty on even for low-mass WIMPs (see Fig. 10). 
For small M^, since {Eyi) oc VqM"^ (from Eqn. 19), 



^^0 



(33) 



so systematic uncertainties on the WIMP mass due to halo uncertainties 
are of order 20%. Ultimately, sufficient measurement of the energy spec- 
trum may allow better determination of (especially if the WIMP mass is 
determined independently from collider data) and even identification of the 
full WIMP velocity distribution, and hence the shape of the dark matter 
halo. The right panel of Fig. 10 shows the differences in the expected spec- 
trum due to different halo models. Although detectors that are sensitive 
to the direction of the WIMP are in the prototype stage (see Section 3.5), a 
high-statistics detection with a detector capable of the determining the re- 
coil direction would allow the detailed determination of the WIMP velocity 
distribution,''^' essentially ushering in an age of WIMP astronomy. 

Finally, better measurements of WIMP mass from colliders may be 
combined with information from direct detection to better constrain the 
WIMP-nucleon cross section (and hence particle-physics parameters). For 
many models, the LHC will constrain the WIMP mass to 10%. However, 
it is difficult to measure WIMP properties well. If the LHC determines the 
WIMP mass, direct detection can determine WIMP-nucleon cross section 
much better than LHC alone. 
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Recoil Energy (keV) Recoil Energy (keV) 

Fig. 10. Left: Expected interaction rate on Ge as a function of recoil energy for different 
WIMP masses (as shown in the legend in units of GeV/c^), for most probable WIMP 
velocity dq = 270km/s (solid) and vq = 170km/s (dash), which span the extremes of 
possible values. For a large-statistics detection of a low-mass WIMP, the 30% uncertainty 
on VQ would result in a similar uncertainty on the WIMP mass, as seen by the similarity 
in spectral shapes e.g. between the lOGeV/c^ WIMP with vq = 270km/s and the 
lOOGeV/c^ WIMP with dq = 170km/s. Right: Expected interaction rate on Ge for 
a 60-GeV/c^WIMP as a function of recoil energy for isothermal, triaxial, and Evans 
halo models with various parameters. For an initial WIMP discovery, uncertainty in the 
WIMP velocity distribution will increase uncertainty on the WIMP mass inferred from 
the detection. Ultimately, these differences in energy spectra may allow inference of the 
correct halo model. 

3. WIMP Direct Detection Experiments 

Direct-detection experiments have already limited the expected WIMP- 
nucleon interaction rate to fewer than 1 event per 10 kilograms of target 
material per day (10 kg-day). With such a small event rate, it is a daunting 
task to search for a WIMP interaction amongst the background interactions 
from cosmic rays and natural radioactivity, which typically number in the 
millions per kg-day (see Ref. 76 for a detailed review of the principle sources 
of background for underground experiments) . 

Because it is not possible to distinguish a single neutron scatter from a 
WIMP scatter if the neutron does not scatter in additional active material, 
neutrons provide a particularly dangerous background for WIMP-search 
experiments. Material with lots of hydrogen, such as polyethylene or clean 
water, acts as shielding for neutrons by reducing the neutrons' energies 
enough that they cannot cause a recoil above threshold (due to simple 
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kinematics, more massive elements do not significantly reduce the energy 
of a a scattering neutron). Neutrons produced by {a,n) reactions (from 
uranium and thorium in rock walls, for example) may be effectively shielded 
in this way since these neutrons start with relatively low energies and have 
high interaction cross sections. For every 13 cm of polyethylene, this low- 
energy neutron background is reduced by an order of magnitude, with a 
thickness of ^40 cm or more needed for current experiments. 

Such shielding is not effective for more energetic neutrons, such as those 
produced by cosmic-ray muons. To reduce this critical background (and oth- 
ers from cosmic rays), all experiments are located underground, with all but 
prototype experiments located deep underground. Table 3 lists the depths 
and locations of the principle underground laboratories for dark matter ex- 
periments. Since denser rock provides a greater effective depth than less 
dense rock, depths are standardly listed in terms of the thickness of water 
{e.g. meters of water equivalent, or mwe) that would provide the same inte- 
grated density as the actual overhead rock. For facilities under mountains, 
usually the mean effective depth is quoted, which inaccurately suggests a 
lower muon flux than in actuality, since shorter pathlengths dominate the 
muon flux. As shown in Fig. 11, depth is effective for reducing backgrounds 
due to the energetic neutrons produced by cosmic-ray muons. In addition, 
most experiments are surrounded, or at least covered, by an active muon 
veto to allow rejection of energetic neutrons if the muon progenitor passes 
close to the experiment. Designs for most future experiments use large in- 
strumented water tanks to provide both shielding for low-energy neutrons 
and identification of fast neutrons or muons that traverse the shield. 



Table 3. Locations, depths, and eflfectivc dopths'^^ of primary underground facilities 
for dark matter experiments.^® Both the 4850-foot (currently the Sanford Lab) and 
planned 7400-foot DUSEL spaces arc listed. 





Depth 


Depth 




Laboratory 


(m) 


(mwe) 


website 


WIPP, AZ, U.S. 


600 


1600 


www.wipp.energy.gov/science/index.htm 


Soudan, MN, U.S. 


710 


2000 


www. hep . umn.edu /soudan 


Canfranc, Spain 


<860 


<2450 


www.unizar.es/lfnae 


Kamioka, Japan 


1000 


2050 


www-sk. icrr.u-tokyo.ac.jp 


Boulby, U.K. 


1100 


2800 


hepw WW .rl.ac.uk/ukdmc/pix /boulby. html 


Gran Sasso, Italy 


1400 


3100 


www .lngs.infn.it 


Modane, France 


1760 


4200 


www-lsm.in2p3.fr 


Sudbury, Canada 


2160 


6000 


www.sno.phy.queensu.ca 


DUSEL, SD, U.S. 


1500 


4300 


www.dusel.org 


(planned) 


2260 


6200 
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Fig. 11. Flux of muons (dashes) and muon-induced neutrons (solid) as functions of 
depth underground, measured in terms of the equivalent thickness of water in km below 
a flat surface that is needed to provide equal shielding. Effective depths of primary un- 
derground facilities for dark matter experiments arc listed (Canfranc is similar to Soudan 
or Kamioka) . Although the neutron background resulting from a given fast neutron flux 
is highly dependent on the experimental setup and materials, the curve of neutron flux 
(still solid) referred to the right-hand axis shows the limit on sensitivity reach due to 
neutron backgrounds for one possible experimental setup.'^^ At depths below about 10 
km w. e., the muon flux is 3 X 10~^^ cm~^ s~^, dominated by neutrino-induced muons. 



Experiments take additional precautions against other sources of back- 
grounds, which otherwise would cause ~ 10'* events keV^^kg^^ day" ^. Low- 
radioactivity copper, which is straightforward to produce, or lead with an 
inner liner of ancient lead (for which radioactive isotopes present at its 
smelting have decayed away) is used to reduce the background from pho- 
tons, typically by 4-5 orders of magnitude. Any air near the detectors is 
purged of radon. Figure 12 shows a typical shielding setup around an exper- 
iment. Materials that surround or constitute the detector must be ultra- 
low-radioactivity, requiring they be screened for possible contamination. 
Residual radioactivity in the detectors or their shielding typically is the 
dominant source of background in experiments, with radioactivity on de- 
tector surfaces (typically from plateout of radon daughters) a particular 
problem. 

In addition to shielding backgrounds, experments reject events that 
are more likely to be due to backgrounds such as photons, electrons, or 
alpha-particles. For example, WIMPs interact so weakly that they never 

interact more than once in a detector, allowing experiments to reject 
multiple-scatter events. Most detectors allow rejection of some multiple- 
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Plastic Detector Ancient Dilution 




Fig. 12. Sliced side view of a typical shielding setup (liere for tiie CDMS experi- 
ment^*). Outermost scintillator paddles act as a veto against events due to muons. 
A thicli polyethylene shield moderates the fiux of neutrons from {a, n) reactions by 5 
orders of magnitude. Lead (grey) reduces the photon backgrounds by 4 orders of mag- 
nitude, with the inner ancient lead liner (dark grey) reducing the background from 
electron bremsstrahlung from ^^''Bi, a daughter of ■^^'^Pb present in modern lead. In- 
ner polyethylene reduces the neutron background from fast neutrons that penetrate 
the outer polyethylene and interact in the lead. Additional polyethylene within the lead 
shield would improve the neutron moderation, but at a significant increase in the amount 
and expense of the lead shielding. Materials that cannot easily be made radiopure (here 
the dilution refrigerator used to cool the cryogenic detectors) must be shielded from the 
detectors. 

scatter events with negligible loss in efficiency to WIMPs, often through 
the use of arrays of detector modules, so that if two separate modules have 
energy depositions, the event must be a multiple-scatter. For liquid nobles, 
multiple scatter may result in pulses separated in time or energy deposited 
in places sufficiently separated in space so as to allow identification by event 
reconstruction. Similarly, WIMPs interact uniformly throughout a detector, 
so it pays to cut interactions near detector surfaces, where more background 
interactions occur. Most experiments use some form of event reconstruction 
to form a "fiducial" volume by rejecting events inferred to occur near the 
detector surface. When comparing detector masses, it is most appropriate 
to consider this fiducial mass. 

Most significantly, WIMPs tend to interact with an atom's nucleus, 
while the dominant radioactive backgrounds (everything except neutrons) 
interact with electrons, so experiments that discriminate between interac- 
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tions causing an electron to recoil and those that cause a nuclear recoil can 
reject virtually the entire radioactive background. There are three ways to 
discriminate between electron recoils and nuclear recoils. Each is based on 
the fact that, for « lOkeV, an electron recoils with v w 0.3c, whereas a 
nucleus recoils with w « 7x lO^^c, depositing its energy much more densely 
over a very short track. Threshold detectors such as COUPP''^ (described 
in more detail in Section 3.1) require a dense energy deposition to trigger 
and therefore are nearly immune to electron recoils, whose deposited en- 
ergy is almost never dense enough to trigger. For other experiments, the 
pulse timing is different for electron recoils than for nuclear recoils. Fi- 
nally, depending on the material, recoil energy may be converted into light, 
ionization, and/or phonons. Experiments that measure two of these forms 
may discriminate against electron-recoil backgrounds because the relative 
amount of energy in the two forms is different for nuclear recoils than for 
electron recoils. 

These different measures of recoil energy and the differing response of 
electron and nuclear recoils may introduce an ambiguity in quoted energy. 
To avoid this ambiguity, most (but not all) experimenters are explicit about 
which signal (or combination of signals) is used to determine an event's 
energy. The unit "keVee" quantifies a measured signal in terms of the energy 
(in keV) of an electron recoil that would generate it, while "keVr" indicates 
the energy of a nuclear recoil that would generate the signal observed. 

The energy scale for keVee is generally easy to establish since photon 
backgrounds (or calibration sources) typically produce mono-energetic fea- 
tures at known energies, although sometimes extrapolations of the scale to 
low energy arc required. The nuclear-recoil energy scale is more difficult, 
due to a lack of such features. In practice, two methods are used. Neu- 
tron scattering experiments allow the nucleus's recoil energy to be inferred 
from the incoming neutron energy and the neutron's measured angle of 
scattering, while simultaneously measuring the signal size. Unfortunately, 
multiple-scatter backgrounds are usually bad enough that significant simu- 
lations are needed to obtain accurate results, especially for the low energies 
that can be measured only with forward scattering. Alternatively, compar- 
ing simulation results to the observed shape of an energy spectrum from a 
neutron source with a broad energy may yield the energy scale. Often, the 
ratio of a signal in keVr to keVee is called the signal's quenching factor, 

QF = E{keYT)/E{keYee). (34) 



(although sometimes the inverse of this quantity is called the quenching 
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factor). If an electron recoil would produce a larger signal than a nuclear 
recoil of the same energy (as is usually the case) , the quenching factor of 
the signal QF > 1. 

There are large fundamental differences in light, ionization, and phonon 
signals. Light signals are the fastest, with ns timing possible, but only ^10 
photons are produced per keV. In order to take advantage of the excellent 
discrimination potential of timing using light signals, efficient light collec- 
tion is critical. Ionization is somewhat better, with ^100 quanta per keV, 
while a whopping ~10,000 phonons are prodiiccd per keV. Experiments that 
detect phonons therefore have fundamentally better energy resolution and 
energy-based discrimination capabilities compared to other experiments. 

Reduction of backgrounds is critical in order to maximize sensitivity 
reach. If backgrounds are kept negligible, the search sensitivity of an ex- 
periment is directly proportional to the target exposure (target mass Mx 
exposure time t). If the expected background is non-negligble but can be 
estimated with negligible uncertainty by some means, it may be statistically 
subtracted (explicitly or implicitly), with the resulting Poisson errors caus- 
ing the sensitivity improvement to be proportional to VMT.^° In practice, 
most dark matter experiments have been background-dominated, without 
means to estimate the backgrounds accurately. In these cases, the experi- 
ments are unable to take full advantage of their target exposure, as increas- 
ing exposure would result in little or no sensitivity improvement due to the 
systematic uncertainties in any background subtraction. The importance 
of systematics also makes pursuing different techniques critical, since these 
technologies tend to have different systematics, thus providing critical cross 
checks for a detection claim. 

The basic techniques include threshold detectors that nucleate if a suffi- 
cient energy deposition occurs (see Section 3.1), ultrapure scintillators (see 
Section 3.4), masses of liquid nobles view by light detectors, with or without 
an electric field to collect ionization (see Section 3.2), solid-state detectors 
cooled to mK temperatures in order to detect phonons, as well as light or 
ionization (see Section 3.3), gaseous detectors for measuring the direction 
and energy density of each recoil (sec Section 3.5), and others that do not 
fit into these categories. There are dozens of WIMP-search experiments in 
progress or development worldwide. Table 4, based on the more extensive 
but slightly dated table in Ref. 81, lists a selection of them. In the sections 
that follow, I describe the various types of WIMP-search experiments, with 
the aim of introducing how the technique works and describing its advan- 
tages and challenges. I have not included discussion of many interesting re- 
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Table 4. Characteristics of selected dark matter experiments,*^ including fiducial mass 
M and whether scintillation light (7), phonons (0), ionization (q), or another form 
of energy is detected, and whether the experiment's primary mission is neutrinoless 
double-beta decay (/3/3). 







Readout 
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Search 


Experiment 


Location 


(t, ?) 


(K) 


(kg) 


Target 


Dates 


NAIAD 


Boulby 


7 


300 


50 


Nal 


2001- 


-2005 


DAMA/Nal 


Gran Sasso 


7 


300 


87 


Nal 


1995- 


-2002 


DAMA/LIBRA 


Gran Sasso 
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300 
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Nal 
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2UUo- 
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Canfranc 
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300 
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Nal 


2000- 


-2005 


ANAIS 
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300 
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2011- 
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Yangyang 
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Yangyang 
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suits that have been made pubhc during the preparation of this manuscript. 
For more detailed reviews of dark matter experiments, see Refs. 81-85. 

3.1. Threshold Detectors 

As mentioned above, one promising technology uses a superheated liquid 
(in bulk^^ or as droplets within a matrix^*'-^^) as a threshold detector. By 
tuning thermodynamic parameters {e.g. temperature and pressure), the 
detector may be made insensitive to the low energy density deposited by 
a minimum-ionizing electron recoil. Only a dense energy deposition, such 
as from a nuclear recoil, will provide enough energy to cause nucleation 
(smaller depositions result in sub-critical bubbles that are squashed to noth- 
ing by their surface pressure). A drawback of these detectors is that they 
provide limited energy information on events. The attraction of these de- 
tectors is that they could allow inexpensive scaling to very large masses 
with a broad range of materials and without need of cryogens or photon 
shielding. Already, the ability of these experiments to use; materials sen- 
sitive to spin-dependent interactions on protons has allowed them to set 
world-leading limits, ^"^'^^'^^ as shown in Fig. 2. 

The primary background for these experiments other than neutrons 
(which plague all experiments) is alpha decays from radioactive contam- 
ination (primarily from radon daughters) in the detectors. The PICASSO 
collaboration has recently shown that these alpha decays are rejectable.*^ 
Since a recoiling alpha particle is lighter and hence has a much longer stop- 
ping distance than a recoiling nucleus, an alpha particle may produce more 
than a single nucleation along its track. Although these individual nucle- 
ations grow together quickly, the amplitude of the high-frequency signal, 
detectable through microphones on the detector, allows discrimination from 
the nuclear recoils of WIMPs. If the rejection proves to be strong enough, 
it will dramatically increase the achievable sensitivity of experiments based 
on the superheated liquid technique. 

One of these experiments, COUPP,^^ uses a large volume of super- 
heated liquid and hence is similar to the classic bubble chambers used in 
high-energy experiments in the 1960s. A principle challenge is keeping the 
detector stable, since any nucleation in the chamber requires recondensing 
the liquid, with a significant deadtime. The COUPP collaboration has met 
the principle challenge of preventing nucleation due to micro-cracks in ves- 
sel walls. They are currently running a 60-kg chamber with higher-purity 
materials and capability of rejecting alpha events, and have shown prelim- 
inary results indicating an impressive reduction in backgrounds, including 
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Fig. 13. Diagram of the processes leading to primary scintillation ("SI") light in a 
liquid noble detector (here Xe), and (if the detector is dual-phase) to secondary ("S2") 
light proportional to the amount of ionization. Recoils dissipate energy as atomic motion, 
excitation, and ionization. Both excitation and ionization result in excited dimers, Xcj, 
in either a longer-lived triplet state or a shorter-lived singlet. Numerical values are for 
Xe; see Table 5 for properties of Ar and Ne. Figure based on Ref. 91. 

a significant reduction in the rate of alplras from impurities in the quartz 
of the vessel that otherwise would limit the livetime of a larger device. 



3.2. Noble Liquid Detectors 

Detectors using noble liquids (and/or gases) also show great promise for 
WIMP detection and have the advantage of relatively easy scaling to large 
masses. Figure 13 shows the basic physics behind these detectors. A recoil in 
liquid Xe, for example, induces both ionization and excitation of Xe atoms 
(in addition to wasting some energy increasing atomic motion). Both exci- 
tation and ionization lead to production of either a singlet or triplet state 
of an excited dimer (Xej). De-excitation of either state produces emission 
of a 175-nm photon that is not absorbed by the noble liquid. The triplet 
state has a longer lifetime than the singlet state (27 ns vs. 3 ns for Xe). The 
dense energy depositions from nuclear recoils result in fewer triplet decays 
and faster recombination, so nuclear recoils have a faster pulse shape than 
electron recoils. As shown in Table 5, the effect is particularly pronounced 
in Ar and Ne, leading to extremely good discrimination in Ar and Ne based 
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Table 5. Properties of noble liquid detectors.®^'' 



Property (unit) 


Xe 


Ar 


Ne 


Atomic Number 


54 


18 


10 


Mean relative atomic mass 


131.3 


40.0 


20.2 


Boiling Point Tb (K) 


165.0 


87.3 


27.1 


Melting Point Tm (K) 


161.4 


83.8 


24.6 


Liquid density at T^, (g cm^'^) 


2.94 


1.40 


1.21 


Volunie fraction in Earth's atmosphere (ppm) 


0.09 


9340 


18.2 


Cost/kg'-^ 


$1000 


$2 


$90 


Scintillation light wavelength (nm) 


175 


128 


78 


Triplet lifetime (ns) 


27 


1600 


15000 


Singlet lifetime (ns) 


3 


7 


<18 


Electron mobility (cm'^ V^^ s^^) 


2200 


400 


low 


Scintillation yield (photons/keV) 


42 


40 


30 



Note: ^ price subject to change; the author does not guarantee any price 
listed. 



on timing alone. Additional discrimination is possible in Xe or Ar based 
on the relative amount of primary scintillation versus ionization (low elec- 
tron mobility makes measuring ionization in Ne impractical). For a given 
energy deposition, nuclear recoils produce less ionization than electron re- 
coils, and much of the ionization of nuclear recoils is "quenched" through 
recombination. 

Experiments with noble liquids must overcome significant challenges. 
Since there are relatively few light quanta (typically 5-10 pe/keV collected), 
maintaining discrimination to low energies is difficult. Stringent purity lev- 
els (~ 10~^ impurities) must be achieved in order to prevent absorption of 
the scintillation light or attachment of drifting electrons. The experiments 
must overcome radioactive backgrounds such as ^^Kr in Xe or especially 
^^Ar in liquid Ar (which produces 1 decay per second per kg of natural 
Ar). 

There are two basic types of noble detectors: single-phase (usually liq- 
uid), which detect only the primary ("S'l") light signal, and dual-phase 
time-projection chambers, which employ a large electric field to drift ion- 
ization electrons upwards out of the liquid and into a region where the 
noble is in its gas phase (see Fig. 14). There, the electrons produce a 
large secondary ("52") light signal by electroluminescence that is propor- 
tional to the amount of ionization. The ratio of this secondary light to 
the primary scintillation {S2/S1) provides additional discrimination with 
typical background leakage 10~^-10~'^ at 50% acceptance. With Xe, dual- 
phase detectors are the only way to get appreciable discrimination against 
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Fig. 14. Left: Sketch of a typical modern single-phase detector geometry. A spherical 
volume of noble liquid is surrounded by an array of close-packed photomultiplier tubes 
(PMTs) to maximize light collection and simplify position reconstruction. A "neck" leads 
up from the top of the detector to allow insertion of calibration sources. Right: Sketch of 
a typical dual-phase detector geometry. Noble gas lies above the liquid with transition 
from liquid to gas occurring at the grid. PMTs on top or bottom of the cylinder allow 
mm reconstruction accuracy in the x-y plane of the detector. Interaction in the sensitive 
volume between the cathode and the grid produces primary scintillation ("5i") light and 
ionization electrons. These electrons drift upwards in the strong electric field produced 
by the anode and cathode, producing secondary ("S2") light proportional to the amount 
of ionization, due to electroluminescence in the gas. As depicted, the 52 signal is much 
larger than the Si signal. The drift time of the electrons causes a delay between the Si 
and S2 signals that allows mm reconstruction of the position in the vertical axis of the 
detector. Interactions occurring in liquid below the cathode produce Si signals but no 
S2 signal, since the direction of the electric field drifts electrons from these interactions 
to the bottom of the detector. 

electron-recoil backgrounds since timing discrimination is weak. However, 
the primary advantage of the dual-phase set-up is that drifting the electrons 
yields mm-accurate position information on the interaction (compared to 
cm-accurate for single-phase detectors), providing much better rejection of 
events due to contaminants external to the detector or on its walls. There 
are costs to these advantages. Drifting electrons is very slow (~ 50 /is) 
compared to collecting scintillation light. Single-phase detectors may stand 
a higher event rate, important for calibration and increasingly important 
as detectors become larger. Since they don't need an electric field, single- 
phase experiments can achieve better light collection, allowing lower energy 
thresholds and/or better pulse-shape discrimination. Their simple design 
allows easiest scaling to very large masses. In contrast, noble liquid exper- 
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Fig. 15. Results of calibration of the DEAP-1 detector with a ^^Na gamma-ray source at 
Queen's University, demonstrating > 10^:1 rejection of electron recoils based on their low 
fraction of prompt light Fprompt- Nuclear recoils have median -Fprompt = 0.82 (dashed 
vertical line labeled "50% acceptance"). Plotted is the probability Pieak that an electron 
recoil would have a fraction of prompt light greater than -Fprompt • 

None of the 1.53 X lO'^ 

electron-recoil events (circles with error bars indicating uncertainties due to counting 
statistics) have Fprompt > 0.64, safely below the value needed to accept 90% of WIMPs 
(dashed vertical line labeled "90% acceptance"). Data show agreement with predictions 
from simulations when a statistical model of photoelectron noise is included (shaded 
swath); the model without photoelectron noise is also shown (dotted curve). The full 
model predicts ~ 10® :1 rejection of electron recoils at 90% acceptance of WIMPs. Cal- 
ibration is underway deep underground at SNOLAB to measure rejection beyond that 
demonstrable at Queen's due to expected neutron backgrounds (solid curve). Figure 
courtesy of C. Jillings. 



iments that measure ionization are inherently not as scalable since large 
detectors would require extremely high voltages and high purity for effi- 
cient collection of electrons; Ar detectors that measure ionization would be 
limited by the high event pile-up due to the relatively slow electron drift 
speed. 

Three single-phase experiments are under construction and should take 
data in 2010 or 2011. XMASS,^* which uses Xe, will take advantage of self- 
shielding to create a low-background 100-kg fiducial volume within 800 kg 
of instrumented Xe. The collaboration appears to be able to achieve suf- 
ficient reduction of backgrounds in the bulk Xe, most notably reduction 
of ^^Kr to 3ppt by distillation.^^ Experiments of the DEAP/CLEAN col- 
laboration take advantage of the large timing difference between nuclear 
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recoils and electron recoils in argon or neon. As shown in Figure 15, the 
DEAP-1 experiment has already demonstrated better than 10'':1 rejection 
of backgrounds,^^ and an analytic model predicts better than billion-to- 
onc rejection with 90% efficiency for WIMPs. Construction of the 100-kg 
fiducial-mass MiniClean dctcctor*^^ and the 1000-kg fiducial-mass DEAP- 
3600 detector^® are both in progress and should start taking data in 2011- 
2013. 

As shown in Figure 1, the XENON-10 experiments^ has the most sensi- 
tive current results among noble-liquid experiments,'*^ including the most 
stringent spin-independent limits of all experiments for WIMP masses be- 
low 40GeV/c^, and, as shown in Fig. 2, the most stringent neutron-spin- 
dependent limits of all experiments.^^ These results are from a "blind" 
analysis, in which data-selection cuts including energy range were all set 
before looking at the parameters of candidate WIMP events in the data. 
Such blind analyses arc critical for dark matter searches when data-selection 
cuts are made in a large dimensional space, otherwise potentially allowing 
bias to lead to fine-tuned cuts that preferentially select or omit the candi- 
date events with little effect on the reported expected background leakage 
or WIMP efficiency. This analysis yields about 50% efficiency for nuclear 
recoils for 58.6 days exposure of a detector of 5.4 kg fiducial mass. The 
expected background based on the expected Gaussian tails of the electron- 
recoil log(S'2/6'i) distribution was T.Oj^^ J events. An additional, unesti- 
mated background resulted from multiple-scatter events where one scatter 
occurs in the sensitive volume and the other in the liquid below the cathode 
(see Fig. 14). Since the latter interaction produces light but no ^2 light, 
these events have reduced log(S'2/S'i) values. In the blindly chosen energy 
region, 2-12keVee (corresponding to 4.5-27 keVR for the assumed nuclear 
quenching factor, but to a higher energy range for more recent measure- 
ments^-^'-^'^''), 10 presumed background events passed cuts (see Fig. 16), with 
studies after unblinding suggesting that most were multiple-scatter events 
with one interaction below the cathode. 

A similar experiment, ZEPLIN-III,^*^^ obtained very similar results^*'^^ 
with a different geometry. While XENON-10 was relatively long and thin, 
with a 20 cm diameter and a 15 cm maximum drift length resulting in a 
drift field of 0.7kVcm~^, ZEPLIN-III is short and squat, with a diam- 
eter of 38.6 cm and a maximum drift length of 3.5 cm, allowing a larger 
drift field of 3.9kVcm^* and better resulting S2/S1 discrimination. ^^'^"^ 
Their analysis was not blind, instead including strong efforts to remove the 
multiple-scattering events that limited XENON-10. A difference in photo- 
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Fig. 16. Plots of events observed with dual-phase Xe detectors. Left: Corrected ion- 
ization yield (which is log{S'2/5'i) rescaled so that photons lie at log(S'2/Si)corr = 2.5) 
and primary scintillation Si energy for XENON-10 events passing blind data-selection 
cuts. 10 WIMP candidate events (circled) lie between 2— 12keVee with ionization yield 
between the median and — 3cr values for nuclear recoils (horizontal lines as function of Si 
bin), consistent with backgrounds from electron-recoils including multiple-scatter events. 
Plot from Ref. 47. Right: Ionization yield log(S'2/5i) and primary scintillation Si en- 
ergy for ZEPLIN-III events passing data-selection cuts. Curves depicting the median 
(central, red curve) and ztltr (outer, blue curves) ionization yield for nuclear recoils are 
shown. 7 WIMP candidate events (large filled circles) lie between 2— 16keVee with ion- 
ization yield between the median and —2cr values expected for nuclear recoils. All are 
near the top edge of the nuclear-recoil candidate region, consistent with expectations 
from background. Plot from Ref. 48. 



tube performance between low-rate WIMP-search running and high-rate 
cahbrations made estimation of expected backgrounds difficult, but the 7 
observed events all appeared near the edge of the S2/S1 cut position and 
so appeared consistent with background, as shown in Fig. 16. Limits for 
XENON-10 arc significantly better at low WIMP masses than those for 
ZEPLIN-III because XENON-10 observed no events between 2-7 keVee, 
and ZEPLIN-III assumed a much more conservative scintillation efficiency 
for nuclear recoils in Xe. 

Two follow-up experiments (50-kg fiducial-mass XENON-100^"^ and 
100-kg LUX^"'^) should have science results in 2010 or 2011, although con- 
struction of LUX has been greatly delayed due to the need to drain the 
flooded Sanford (future DUSEL) site. If backgrounds are well below the 
rate seen in XENON-10, sensitivity ~10~^ pb may be achieved. Prototypes 
for dual-phase Ar detectors WArP"*^ and ArDM^""* have also been con- 
structed. Use of argon with low ^^^Ar content^"^ would allow operation of 
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massive detectors without event pile-\ip. Running the 3.2 kg fiducial WArP 
prototype resulted in no candidate events above 42keV in a 96.5 kg day 
exposure. The WArP collaboration is building a 140-kg experiment that 
could have results in 2010 or 2011. The ArDM collaboration is constructing 
a 1-ton detector at CERN and will consider a deep underground opera- 
tion after its successful commissioning. ^'^^ Relatively poor light collection 
appears likely to limit the detector to a 30keV energy threshold without 
reach beyond that of current experiments. Plans exist for rapid increases 
in mass to ton or multi-ton liquid noble detectors. 



3.3. Cryogenic Detectors 

Ultrapure Ge semiconductor detectors operated at the temperature of liq- 
uid nitrogen, 77 K, were used for the first searches for dark matter particles 
in the 1980s. ^^^'^^^ Because these detectors measured only the ionization 
of energy depositions, they had no discrimination between electron-recoil 
backgrounds and nuclear recoils. Although cryogenic detectors with strong 
background discrimination are now used more widely, ionization-only detec- 
tors with improved designs are operated today because of the very low en- 
ergy thresholds they may achieve. The TEXONO collaboration has achieved 
an energy threshold of 220cVce in four 5g detectors, ^'^"^ while the CoGeNT 
collaboration achieved nearly as low a threshold for a run of 8.4 kg days 
exposure. 

The key advantage of most cryogenic detectors is the ability to measure 
the energy deposited as phonons, vibrations of the crystal lattice. The ini- 
tial phonons produced are not at equilibrium. Some detectors collect these 
"athermal" phonons, which contain information on the location and type 
of recoil that occurred. On timescales of ms, essentially all other forms 
of energy depositions are converted to heat and the phonons thermalize, 
resulting in a temperature increase of the detector AT = E/C ^ l/^K, 
where E is the energy deposited and C is the detector's heat capacity. 
Detectors insensitive to athermal phonons measure this temperature in- 
crease with sensitive thermistors, such as neutron-transmutation-doped Ge 
or transition-edge sensors (see e.g. Ref. 115). 

The EDELWEISS experiment uses thermistors attached to Ge crystals 
at cryogenic temperatures (20 mK) to measure phonons, in addition to mea- 
suring ionization using a small applied electric field. Backgrounds have been 
dominated by low-energy electrons interacting near the surfaces of the de- 
tcctors^^^'^^^ because such interactions result in incomplete collection of the 
ionization charges, mimicking nuclear recoils. The problem occurs because 
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Fig. 17. Sketches of the EDELWEISS interdigitated and CRESST cryogenic detectors. 
Left: Side view of a cylindrical interdigitated detector. Both top and bottom sides have 
two sets of electrodes, with each commonly labeled electrode (a,b, c, or d) connected 
in parallel. The voltage biases chosen {e.g. Vk = 2V, Vb = —IV, Vc = — 2V, = 
1 V) produce an axial field in the bulk of the detector, perpendicular to the depicted 
equipotential lines (with voltages labeled by the vertical axis), while the field close to 
the surface links two adjacent electrodes and is therefore approximately parallel to the 
surface. A bulk event (leftmost in figure) is thus identified by the collection of electrons 
and holes on the "fiducial" electrodes a and c. An event near enough to the surface 
for charge collection to be incomplete (middle or rightmost in figure) would result in 
charge collection in one of the two "veto" electrodes b or d, allowing its rejection. Figure 
based on Refs. 119,120. Right: Sketch of CRESST phonon and scintillation detector. 
Phonons produced in the CaW04 target are detected by the attached transition-edge 
sensor, which is connected by a thermal impedance to a constant-temperature cryostat, 
as in other phonon detectors. Light emitted by the CaW04 target is detected by an 
adjacent silicon-on-sapphire detector, also with an attached TES connected by a thermal 
impedance to the cryostat. Figure from Ref. 121. 



the electrons and holes generated by an interaction are sufficiently ener- 
getic to diffuse against the applied electric field into the nearby electrode, 
causing a fraction of the event ionization to be "lost" if the interaction is 
near the surface. 

The new EDELWEISS "interdigitated" detectors^^" promise to provide 
excellent discrimination against this dominant background. As shown in 
Fig. 17, each side of the detector is interleaved with oppositely charged 
electrodes. Events in the bulk of the detector result in the collection of 
electrons and holes only on the "fiducial" electrodes, whereas events near 
the surface result in some charge collection in the veto electrodes. Overall, 
calibration indicates rejection of surface events has only 10~^ inefficiency, 
sufficient for 80 kg years exposure at achieved backgrounds. Recent run- 
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ning of an array of ten 400 g detectors for an effective exposure of 144 kg d 
resulted in one nuclear-recoil candidate above 20 keV, leading to the com- 
petitive limits shown in Fig. 1.^^ Continued running through Spring 2010 
is planned, at which point additional detectors with larger fiducial masses 
will be commissioned and run. 

The CRESST experiment also uses cryogenic detectors, but measures 
phonons and scintillation light, which provides as good rejection of surface 
events as it does rejection of photons in the bulk since light production is not 
reduced for events near the crystal surface. The main CRESST detectors 
are scintillating (300 g) CaW04 crystals, with the W providing the heavi- 
est nucleus of any dark matter search. As shown in Fig. 17, a transition- 
edge sensor deposited on the detector surface measures the phonon energy. 
Adjacent to each crystal is a small (2g) silicon-on-sapphire light detec- 
tor equipped with another transition-edge sensor. The energy resolution of 
the phonon measurement is excellent, better than 1 keV over the full en- 
ergy spectrum. The energy resolution of the light measurement is worse, 
~10keVee, due to the small fraction of energy converted to light. Nuclear 
recoils produce far less light than electron recoils, especially for interactions 
with the heavy W nuclei (QF « 40). This difference provides excellent dis- 
crimination between electron recoils and nuclear recoils, but also means 
that only high-energy WIMP events would produce any detectable light at 
all, resulting either in a high energy threshold, or in potential susceptibility 
that something causing phonons but no light {e.g. crystal relaxation) may 
mimic a WIMP signal. 

Following major upgrades to the experiment's shielding, detector sup- 
port, and electronics, the CRESST collaboration ran detectors in 2007 in a 
commissioning phase for CRESST-II, yielding a total exposure of 30.6 kg- 
days of W exposure from two detector modules. Three events consistent 
with W recoils were found in the energy region 10-40 keV, with the high 
end of the energy range determined by the cut-off in the WIMP spectrum 
due to the form factor for W (see Fig. 7). As shown in Fig. 1, resulting 
limits for this small exposure are only about 20 x less constraining than 
the world's best. The cause of these events, and of additional events with 
energies > 40 keV, is unknown and under investigation. 

For the longer timescale, the EDELWEISS and CRESST collaborations 
and others have formed a new collaboration, EURECA,^^^ dedicated to a 
cryogenic experiment at or near a ton of detector mass. Plans are for a 
100-kg experiment in 2012 and a ton scale experiment in 2018. 

The CDMS II experiment uses cryogenic Ge or Si ionization-and-phonon 
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Fig. 18. Sketch of a 10-/jm-long CDMS athcrmal phonon sensor on the surface of a 
much larger Ge or Si target. Phonons produced in the Ge or Si target break Cooper 
pairs in the superconducting Al, producing quasiparticles, which diffuse into the overlap 
region, becoming trapped in the W transition-edge sensor. 

detectors that discriminate against the otherwise dominant surface electron 
recoils by collecting phonons before they thermalize.''® Each quadrant of the 
detector's top surface includes thousands of phonon sensors connected in 
parallel. As shown in Fig. 18, athermal phonons produced by the interac- 
tion propagate to the detector surface, where most of them are absorbed 
in superconducting aluminum pads. Quasiparticles generated in the alu- 
minum by the phonons breaking Cooper pairs diffuse in ~10 /is through the 
aluminum to a tungsten transition-edge sensor. The aluminum fins allow a 
large phonon collection area while keeping the transition-edge sensors small, 
so that even a small amount of energy produces a large temperature change. 
Each tungsten sensor is kept in the middle of its sharp superconducting- 
to-normal transition, so a small increase in temperature greatly increases 
its resistance. As a result these detectors have excellent energy resolution. 
Comparison of phonon-pulse arrival times in the four independent chan- 
nels allows localization of the interaction position in the xy-plane of the 
detector. 

The ionization yield allows near-perfect separation of nuclear recoils 
from bulk electron recoils, and the shape, timing, and energy partition of 
the phonon pulses allow rejection of events occurring near the detector 
surfaces. This rejection works because the athermal phonons from electron 
recoils are faster than those from nuclear recoils, particularly if the electron 
recoils occur near a detector surface. Accepting only events with both slow 
phonon pulses and low ionization yield rejects over 99.5% of the surface 
events while keeping over half of the nuclear-recoil events in an analysis 
tuned to maximize the discovery potential of the search. 
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Fig. 19. Normalized ionization yield (number of standard deviations from the mean 
nuclear-recoil yield) versus normalized timing parameter (timing relative to acceptance 
region) for WIMP-search events (+) and neutron calibration events (grey dots) from the 
final exposure of the CDMS II experiment. Only events between 10— lOOkeV, consistent 
with all signal criteria (excluding yield and timing) are shown. Two WIMP candidate 
events lie within the solid box indicating the signal region. Figure courtesy Zeesh Ahmed. 



Data taken with ^4 kg of Ge detectors in 2007-2008 resulted in a total 
exposure to WIMPs of 612kg-days. An analysis in which data-selection cuts 
were set blind, based on events from calibration sources or other events 
that could not be from WIMPs, resulted in an expected background of 
0.8 ± 0.1(stat)±0.2(syst) misclassified surface electron recoils and ^ 0.1 
events from neutrons. Two WIMP candidates, at recoil energies of 12.3 keV 
and 15.5 keV, were observed (see Fig. 19). Because the probability to have 
observed two or more background events in this exposure is 23%, these 
results are not significant evidence for WIMP interactions.'*'* Combined 
with the results of previous exposures, resulting spin-independent limits are 
the world's most constraining for WIMP masses > 40GeV/c^, as shown in 
Fig. 1. 
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A follow-up experiment, SuperCDMS,^^'*'"'^^''' is being commissioned us- 
ing the same infrastructure but improved, thicker (inch-thick instead of cm- 
thick) detectors with reduced exposure to radon daughters, and hence lower 
surface backgrounds. The larger detectors result in ~ 4x the total mass in 
Ge, with sensitivity reach to 4 x 10~^ pb, at which point the neutron back- 
grounds of the current setup at Soudan likely would begin to dominate. An 
experiment with 100 kg at a deeper site, SNOLAB, is planned. Research and 
development on the manufacturability of the 1.5 ton GEODM project, '^^^ is 
in progress. It may be possible to make detectors 20 x larger than the CDMS 
detectors by exploiting inexpensive, dislocation-free Ge crystals, which are 
unusable when run at liquid nitrogen temperatures but appear to work at 
50 mK. Multiplexing would simplify the readout electronics and reduce the 
heatload on the cryogenic systems. 

3.4. DAMA 

The only claimed detection of WIMPs is from the DAMA experi- 
ments. ^^''"^■'^ DAMA does not distinguish between WIMP signal and back- 
ground on an event-by-event basis. Instead, the presence of WIMPs may 
be inferred from the annual modulation in the rate of the lowest-energy 
single-scatter interactions, assuming that the backgrounds do not modu- 
late significantly. The DAMA/Nal apparatus, consisting of 87.3 kg of Nal 
scintillator crystals, was run from 1996-2002. The current, DAMA/LIBRA 
apparatus^^*^ consists of an array of 25 Nal scintillator crystals, each 9.7 kg. 
One is not operational, resulting in a total target mass of 232.8 kg. Each 
crystal is viewed by 2 phototubes through suprasil-B lightguides, with ex- 
cellent achieved light collection (5.5-7.5 photoelectrons/keVee). Extreme 
efforts have been taken to avoid contamination, including etching of parts 
followed by installation within a high-purity nitrogen atmosphere using 
breathing apparatus. Photon and neutron shielding is fairly standard, al- 
though it may be notable that the detectors are separated by a fair mass 
of copper, and there is no surrounding scintillator veto. 

As shown in Fig. 20, the annual modulation in DAMA's event rate is 
compelling. The combined 0.82 ton-years of exposure result in a best-fit 
amplitude A = 0.0131 ± 0.0016 keVee"^kg-May"\ period T = 0.998 ± 
0.003 years, and phase to = 144 ±8 days (consistent with to = 152.5 days as 
expected for a standard halo), with a significance of 8.2cr. As expected for 
standard WIMP elastic scattering discussed in Section 2, the modulation 
amplitude (sec Fig. 21) is significant only at low energies (< 5keVee). The 
shape of the energy spectrum of the modulation is not well constrained at 
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Fig. 20. Sum of residuals of the single-hit scintillation events in the 2-6 keVec energy in- 
terval, after subtracting time-averaged rates in each energy bin in each detector, as a func- 
tion of days since January 1, 1996, for the DAMA/Nal and DAMA/LIBRA experiments. 
The experimental errors are vertical bars, and the associated time bin widths are horizon- 
tal bars. The superimposed curve represents the cosiimsoidal function Acosu!{t—to) with 
modulation amplitude A = (0.0129±0.0016) keVee~^kg~'^ day" ^obtained by best fit over 
the whole data while constraining the period T = 27r/uj = 1 yr and phase to = 152.5 day 
(June 2"'^). The dashed vertical lines correspond to the maximum of the signal (June 
2"''), while the dotted vertical lines correspond to the minimum. Figure from Ref. 131. 
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Fig. 21. Energy distribution of the binned modulation amplitude Sm for the total 
exposure of DAMA/Nal and DAMA/LIBRA. A clear modulation is present in the lowest 
energy region, while Sm values compatible with zero are present just above. Figure from 
Ref. 131. 



low energies; since the lowest-energy bin has a somewhat lower rate than the 
next higher bins, a spectrum that is monotonically falling with increasing 
energy, or one that rises and then dips (such as a monoenergetic line near 
3 keVee) each provides an adequate fit to the data.^^^ Finally, no modulation 
is seen in the rate of low-energy multiple-scatter events, providing evidence 
that the signal is not a simple modulation of the background. 

If interpreted as a standard WIMP interaction on iodine, the results 



January 28, 2011 1:7 WSPC - Proceedings Trim Size: 9in x 6in SchneeTASIPDF 



40 




WIMP Mass (GeV/c^) WIMP Mass (GeV/c^) 

Fig. 22. Comparison of DAM A (3<t) allowed regions as interpreted in Ref. 69 to pub- 
lished (90% C.L.) upper limits on the spin-independent (left) and proton-spin-dependent 
(right) WIMP-nucleon cross section as functions of WIMP mass under the standard as- 
sumptions*^ about the WIMP distribution in the Galaxy described in Section 2. For 

DAMA, interactions of WIMPs with mass > 25GeV/c^arc dominated by interactions 
on 1, while those < 25GcV/c^are dominated by interactions on Na. Upper limits of 
CDMS'"' (solid), XENON-IO"^ (dashes), and KIMS"''' (circles) arc clearly incompati- 
ble with spin-independent interactions on 1, while limits from PICASSO^* (6-pointcd 
stars), COUPP'''' (5-pointcd stars), and KIMS (circles) are clearly incompatible with 
spin-dependent interactions on I. Parts of the regions including the effects of ion chan- 
neling as presented in Ref. 132 (light shade) are allowed at low masses. Almost the 
entire region assuming insignificant channeling (dark shade) for both spin-independent 
and spin-dependent interactions, even at low masses, is excluded by combinations of 
limits of CDMS on Si"3 (dash-dots), CoGeNT"* (dots), and PIGASSO. 

are clearly inconsistent with limits from other experiments, as shown 
in Fig. 22 for spin-independent and proton-spin-dependent interactions 
(neutron-spin-dependent interactions are excluded even more strongly, since 
Na and I isotopes are odd-p). The results are in conflict with other exper- 
iments if they arc due to standard low-mass WIMP scattering on soditim 
instead of iodine, whether the dominant interaction is spin-independent 
or spin-dependent, if the standard halo model is correct. Since alternate 
halo models may produce a larger modulation for the same WIMP rate,^^'' 
a non-standard halo may improve consistency between DAMA and other 
experiments. 

Figure 22 also shows DAMA allowed regions if ion channeling is signif- 
icant in Nal, as presented in Ref. 132. Ion channeling is an effect observed 
experimentally for nuclei sent into lattice from outside. If the nucleus trav- 
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els down the channel between lattice sites, it transfers more of its energy to 
electrons rather than to other nuclei, producing much more light than usual. 
If WIMP-induced recoiling nuclei are channeled, the sensitivity of DAMA 
would be increased especially for low-mass WIMPS since the channeling 
would allow a low-energy recoil to produce enough light to be detectable. 
However, channeling has not been observed, nor is it expected, for nuclei 
that start on a lattice site, as would be the case for nuclei recoiling from 
a WIMP interaction. In a perfect lattice, no nucleus would be channeled 
by the rule of reversibility. Accounting for thermal vibrations, Bozorgnia, 
Gelmini, and Gondolo find channeling's effect on the DAMA energy spec- 
trum is < 1%.^^^ 

Many non-standard interactions have been proposed (e.g.^^^^^^^) that 
may explain the DAMA signal. Since these are described in Neil Weiner's 
contribution to these proceedings,^ they will not be discussed here. Instead, 
it is worthwhile considering in detail whether the annual modulation sig- 
nal could be caused by something other than dark matter, in order to help 
evaluate how seriously alternate WIMP models should be considered. Could 
some background be causing the annual modulation? Such a background 
would have to fulfill the annual modulation characteristics of a WIMP. Its 
rate would have to vary cosinusoidally over the course of the year, with 
period T = lyear and phase to ~ 152.5 days. The modulation should ap- 
pear only in the lowest-energy, single-detector hits, and should produce a 
consistent amplitude, A < 7%, between the Nal/LIBRA experiments and 
between different detectors in the experiments. 

Although these requirements are extensive, they are not as strenuous at 
they might at first seem. The date of expected maximum signal, June 2, 
corresponds roughly (but not exactly) to summer, and there are of course 
many systematic differences between summer and winter. The lowest-energy 
events are the very ones most likely to be affected by a systematic effect, 
and the requirement that the multiple-scatter events not show modulation 
is not a very strong test, since there are not many multiple-scatter events. 
The requirement that the amplitude A < 7% is barely restrictive at all. 

Because the overall, unmodulated background was designed to be lower 
in DAMA/LIBRA than in DAMA/Nal, the consistency check between the 
two setups is interesting. As shown in Fig. 23, the continuum background 
>4keVee is ~2x lower in the LIBRA setup. Below about 3.5keVee, the 
backgrounds are surprisingly similar, with a slightly higher (statistically 
significant) rate in DAMA/LIBRA than in DAMA/Nal in the range 2- 
2.5keVee. Any explanation of the modulation in terms of backgrounds 
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Fig. 23. Overall singlc-liit rate in the DAMA/Nal (higher curve at high energies) and 
DAMA/LIBRA (lower curve at high energies) experiments as a function of the electron- 
equivalent energy. Between 2-3keVee the spectra cross. The small peak (O.Tperkgd 
spread ax;ross l.SkeVee) evident in the DAMA/LIBRA spectrum near SkeVee is pre- 
sumably due to ''"K. Figure based on those in Refs. 129,131. 



would need to be consistent with this behavior of the overah rate. 

As can be seen by eye from Fig. 20, the amplitude of the annual mod- 
ulation is somewhat smaller for the DAMA/LIBRA setup than for the 
DAMA/Nal setup (although the larger uncertainties on the DAMA/Nal 
data make the difference appear larger than it is). Table 6, from Ref. 130, 
compares the fits to the annual modulation with the function A cos — to) 
for the two individual set-ups and together. The worst consistency is for 
the full 2-6 keVee region. DAMA/Nal shows a modulation amplitude nearly 
twice as large as that for DAMA/LIBRA, A = 0.0200 ± 0.0032 in compari- 
son to A = 0.0107 ± 0.0019. This difference is 0.0093 ± 0.0037, about 2.5a, 
which happens by chance 1.2% of the time. While such a probability is 
certainly not so low as to indicate that the two experiments are incon- 
sistent, it does not provide especially strong constraints on the possibility 
that some background may be causing the modulation. In particular, the 
fact that the Nal modulation is larger primarily in the same region of the 
energy spectrum (> 5keVee, and to a lesser extent > 4keVee) for which 
its background is larger makes it possible that the reduction in background 
between DAMA/Nal and DAMA/LIBRA is what has caused the reduction 
in modulation amplitude. 

The DAMA collaboration have taken strong steps to assure and to check 
that no systematic effect could be causing the modulation. No suggested 
background process appears likely to yield modulation as large as the ob- 
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Table 6. Results of fits to AcosLj{t — to) with the residual rates of the single-hit 
scintillation events, collected by DAMA/Nal, by DAMA/LIBRA and by the two ex- 
periments together in the 2- 4, 2-5, and 2-6 keVce energy intervals.-'-^-'- The last column 
shows the C.L. obtained from the fitted modulation amplitudes. 

A (keVee-ikg-iday-i) T = ^ (yr) to (day) C.L. 

DAMA/Nal 

2-4 keVec 0.0252 ± 0.0050 1.01 ± 0.02 125 ± 30 5.0a 

2-5 keVee 0.0215 ± 0.0039 1.01 ± 0.02 140 ± 30 5.5cr 

2-6 keVee 0.0200 ± 0.0032 1.00 ± 0.01 140 ± 22 6.3a 

DAMA/LIBRA 

2-4 keVee 0.0213 ± 0.0032 0.997 ± 0.002 139 ± 10 e.7a 

2-5 keVee 0.0165 ± 0.0024 0.998 ± 0.002 143 ± 9 6.9ct 

2-6 keVee 0.0107 ± 0.0019 0.998 ± 0.003 144 ± 11 5.6ct 

DAMA Combined 

2-4 keVee 0.0223 ± 0.0027 0.996 ± 0.002 138 ± 7 8.3ct 

2-5 keVee 0.0178 ± 0.0020 0.998 ± 0.002 145 ± 7 8.9ct 

2-6 keVee 0.0131 ± 0.0016 0.998 ± 0.003 144 ± 8 8.2ct 



served signal. Although no possibility appears plausible, the most 
likely ones (not discussed in Ref. 140) are from a possible modulation in the 
rate of 3.2 kcV x-rays from potassium, in the rate of some muon-induced 
x-ray, or in the effectiveness of the phototube noise cut. Rcf. 81 describes 
additional possible systematic effects that also do not appear likely to be 
the cause of the modulation. 

Potassium-40 is a common cause of background events in dark matter 
experiments due to its high abundance and 10^-year halflife. About 10% 
of the time, decays by electron capture to an excited state of *°Ar, 
which subsequently decays by emission of a 1461 keV gamma ray. Since 
the daughter Ar atom is still missing a K-shell electron, 3.2 keV emerges 
when an L-shell electron drops down to fill the vacancy. Since low-energy 
X-rays and electrons have short interaction lengths, this low-energy de- 
position essentially always appears in the parent detector where the '"'K 
decays, but the 1461 keV gamma ray has a longer interaction length and 
so may deposit some or all of its energy in the parent detector, other de- 
tectors, or passive material such as copper. The DAMA collaboration have 
counted events with the 1461 keV gamma ray in one detector and 3.2 keV 
in another, and compared the rate to simulations of this process in their de- 
tectors to infer the natural K contamination in their detectors is ISppb.^*^ 
This rate is about right ^'^''^ to be responsible for the small peak in the overall 
DAMA/LIBRA event spectrum (0.45keVee"^kg"May"^over 1.5 keVee) at 
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about SkeVee (see Fig. 23). A 4% modulation in the rate of the accepted 
single-scatter 3.2keVee events could cause DAMA's observed modulation. 
This size modulation appears too large to be able to be caused by a modu- 
lation of the trigger efficiency and analysis cuts,^''" but given the similarity 
of the '"^K background to the modulation signal it would be prudent for 
DAMA to check the stability of triggering, analysis cuts, and overall rate 
of the 3.2keVee events that are coincident with 1461 keV gamma rays. It 
is also barely imaginable that changing detector response or actual move- 
ment of the K (due to changing humidity, for example) could result in a 
modulation of the single-scatter rate if such a change alters the probability 
that the 1461 keV gamma rays deposit any energy in the parent detector. 
Therefore, it would also bo useful for the DAMA collaboration to check 
whether reasonable changes in their simulation could result in a different 
rate of 3.2 keV single-hit events. 

The known physical side reaction whose modulation best matches that 
of DAMA's signal is the annual modulation of the muon rate due to sea- 
sonal changes in the atmospheric temperature. An increase in temperature 
of the stratosphere decreases its density, reducing the chance of meson in- 
teractions, resulting in a larger fraction of mesons decaying into muons. 
From 1991-1994, MACRO measured the annual modulation of the muon 
rate at Gran Sasso,^*^ finding the phase nearly matches that of the DAMA 
signal. The maximum occurs in mid- June, a few weeks later than the best- 
fit phase of DAMA's signal, perhaps consistent within uncertainties. The 
muon modulation also appears less regular than that of DAMA's signal. 
DAMA claims that the modulation is too small for muon-induced neutron 
interactions to be a significant source of background modulation. An al- 
ternate possibility that has not yet been tested is that muon spallation in 
Nal detectors may create a metastable isotope that decays with emission 
of 3keV.^'''' The lifetime of the metastable state would have to be > 500 /xs 
to avoid DAMA's trigger holdoff time. A beam test of Nal could see if such 
a line exists. 

Finally, it is imaginable that an instrumental effect could be modulating 

annually. For example, DAMA analysis includes a cut on the sharpness of 
a pulse that is used to remove PMT noise events, which otherwise would 
pollute the lowest energy bins. It is not possible to know what the tail of the 
PMT noise distribution looks like, since measurement always has particle- 
induced backgrounds, so it is not clear that the cut removes all PMT noise 
events. It is possible that a small number of noise events passes this cut, 
and the number that passes this cut modulates. Although this possibility 
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is poorly motivated, dark current in PMTs depends exponentially on tem- 
perature, and it is possible that DAMA's temperature stability is not good 
enough to prevent a modulation of this potential effect. 

The DAMA collaboration have made some checks of effects related to a 
possible modulation of noise events. They have determined that the total 
number of noise events does not modulate (but they cannot check whether 
the number passing the cut modulates), and they have checked (using a 
high-rate calibration source) that the efficiency of the cut for particle in- 
teractions does not modulate. Studies of the effect of the cut position on 
modulation amplitude could help rule out a possible modulation of noise 
events passing the cut. If different cut positions all showed no modula- 
tion in the particle efficiency (using the calibration source) and affected 
the WIMP-search modulation amplitude in a way consistent with this eflS- 
ciency for particles, it would be strong evidence that the modulation is due 
to particle origin. 

Other experiments may be able to test DAMA's results more directly 
soon, using similar detectors. The KIMS experiment, situated in Yangyang 
Lab, Korea, has already yielded limits^^ from running four 8.7 kg Csl crys- 
tals for a total exposure of 3409 kg d. Twelve detectors, for a total mass of 
104.4 kg, are now installed. The ANAIS collaboration performed tests on a 
Nal crystal to determine its background and test its energy threshold.^*'' 
While a 2keVee energy threshold was achieved, the "^"K background level 
in the Nal crystals was found to be too high to allow a test of DAMA's 
results. Work on producing cleaner crystals is in progress. 

3.5. Directional Detectors 

Although inherently difficult, detection of the direction of scattering 
WIMPs would provide significant additional information. First of all, such 
a detection would provide a clear signature of the non-terrestrial nature of 
the signal, since the motion of the Earth through the Galaxy results in a 

predominant WIMP "wind" in the opposite direction, from the constella- 
tion Cygnus. As shown in Fig. 24, a detector would see this wind vector 
oscillate over the course of a day. Over the course of a year, this direction- 
ality would repeat each sidereal day, hence rapidly going out of phase with 
the solar day. With a detector able to distinguish the full direction of re- 
coils {i.e. both the axis and which ends of the track are the head and tail), 
~ 10 events would be suflacient to determine that the WIMP signal is not 
isotropic. As mentioned above in Section 2.4, directional detection 
would provide excellent information on the WIMP velocity distribution. '''^'''^ 
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Fig. 24. Approximate spin axis of the Earth relative to the WIMP "wind" due to the 
motion of the Earth through the Galaxy. A detector at latitude +60° would see most 
scattering WIMPs come predominantly from the lab's ceiling at one time during the 
day and from the northern side wall twelve hours later, with the pattern repeating each 
sidereal day. Figure based on that in Ref. 148. 



Unfortunately, recoil track lengths are too short in solid or liquid de- 
tectors (~ 0.1 /^m in a condensed matter target) for direct measurement. 
While the possibility of using a detector with non-uniform directional re- 
sponse has been considered {e.g. Ref. 148 and references therein), gaseous 
detectors appear the most promising, despite the difhculty in making them 
massive enough to detect a WIMP signal. In order to make the recoil track 
length long enough to be imaged (~ 1 mm), traditional gas time-projection 
chambers may be used if operated at low pressure (< 100 Torr). As with the 
liquid noble time-projection chambers, position in the x — y plane is deter- 
mined by the position of charges collected (by crossed planes of wires^^^ , 
micropixels,^'^^ micromegas,^^'^ or via detection of electroluminescence^^^), 
while position along the drift direction is determined by drift time. However, 
since no prompt signal is detected, the interaction time is unknown, and 
position information along the drift direction is only relative, not absolute. 

Directional detectors potentially provide the best observables of any 
dark matter experiments. The total charge indicates the energy of the re- 
coil. The tracklength of electron recoils is much larger than that of nuclear 
recoils, so comparison of the observed tracklength to the energy allows 
excellent rejection of electron-recoil backgrounds. The track itself easily in- 
dicates the axis of the recoil. Measurement of the amount of charge released 
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along the track can allow inference of the sense of the direction of the track 
{i.e. can allow the head of the track to be distinguished from its tail). For 
nuclear recoils, energies of interest are below the Bragg peak, so the energy 
deposition is larger at first, then decreases as the recoiling nucleus loses en- 
ergy. Unfortimatcly, increased scattering near the end of the track, as well 
as straggling in the ion drift, can make correct identification of the head 
and tail of a track difBcult.^^^ 

The four directional detectors listed in Table 4 are all in research and 
development phase, as their low masses do not yet permit competitive sen- 
sitivities. Chosen gases typically maximize sensitivity to spin- dependent in- 
teractions. The NEWAGE collaboration has run a small, 10-g test cell with 
micropixel readout underground, achieving the best limits so far of any di- 
rectional detector, albeit still 5 orders of magnitude less sensitive than 
the world's best proton-spin-dependent limits. The MIMAC collaboration 
has run test chambers with micromegas;-'^''''^ much effort has been concen- 
trated on low-energy measurements of the quenching factors for gases of 
interest. Two of the experiments, DRIFT^^^ and DMTPC,^^* have demon- 
strated head-tail discrimination of the track. The DMTCP measurement 
was at somewhat higher energies and tracklengths than needed for a dark 
matter search (> 200 keV),^^^'^^^ while the DRIFT measurement shows 
some discrimination down to 50keV.^^^ 

It remains to be seen if the WIMP-nucleon cross section is large enough 
that one can build gas detectors with enough mass to detect the signal. 
A 1-ton gas detector, with spin-independent sensitivity reach 10~^pb, 
would be about the size of the largest accelerator-beam detectors. To ob- 
tain sensitivity to ^ 10~^°pb would require a detector ^ 5x the size of 
SuperKamiokande, smaller than proposed next-generation neutrino and 
proton-decay experiments. Identifying a technology to achieve adequate 
spatial resolution over very large areas at reasonable cost will be a chal- 
lenge. The benefits of a directional detection certainly warrant continued 
development to attempt to reach this lofty goal. 

4. Conclusions and Prospects 

Currently, CDMS II and XENON-10 are the most sensitive WIMP-search 
experiments, with limits on the spin-independent WIMP-nucleon cross sec- 
tions below 10~'^pb, probing the MSSM region. Threshold detectors and 
KIMS, a room-temperature scintillator, have the best direct sensitivity 

to spin-dependent interactions on protons, with indirect detection exper- 
iments SuperKamiokande and IceCube more sensitive for some models. 
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Fig. 25. Rate of coherent neutrino interactions above a given threshold energy for Ge 
(left) and Xe (right). For both the diffuse supernova (DSNB, large black dashes) and 
atmospheric (ATM, small green dashes) event rates, the sum of all contributing neutrino 
flavors is shown. The significant background from *B neutrinos (red solid curve) will 
prevent experiments from using threshold energies < 3keV for Xe, < 5keV for Ge, 
< lOkeV for Ar, or < 20keV for F. Figures from Ref. 160. 

Truly, a variety of technologies has been successful and continue to show 
promise for the future. Significant improvements by the end of 2010 are 
likely from XENON-100, EDELWEISS, and COUPP, with many experi- 
ments (LUX, WArP, XMASS, SuperCDMS, and MiniCLEAN, in addition 
to XENON-100 and EDELWEISS) likely to surpass 10"^ pb sensitivity for 
spin-independent interactions by 2012. The next five years should see spin- 
independent sensitivity improve to 10~^pb, or possibly even 10^^" pb if 
promising technologies can achieve the necessary background rejection. In 
particular, it is unclear what is the best technology for the multiton scale 
needed to achieve 10^^^ pb sensitivity. 

The ultimate unrejectable, unshieldable background for future WIMP- 
search experiments will be coherent neutrino-nucleus interactions (see 
Fig. 25).^^" A large rate of interactions of solar *B neutrinos will prevent 
ton-scale experiments from using threshold energies < 5 keV for Ge, < 3 keV 
for Xe, < lOkeV for Ar or Si, and < 15keV for Ne or F (although sub- 
traction of the expected energy spectrum would be possible, and rejection 
would be possible for directional detectors. ^®^) For yet larger detectors, at- 
mospheric neutrinos will provide background events over the entire energy 
range of the WIMP search, likely limiting the sensitivity reach of ultimate 
dark matter experiments to no better than ~ 10~^^ pb. 
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A direct detection of WIMPs would warrant follow-up in order to learn 
as much as possible about them. Measurement of the WIMP recoil spec- 
trum with good statistics would constrain the WIMP's mass, potentially 
demonstrating that a particle produced at accelerators indeed comprises 
the dark matter in the galaxy. Such large statistics (provided experimental 
operation is sufficiently stable) could also take advantage of the expected 
annual modulation of the WIMP signal to confirm the extraterrestrial ori- 
gin of the WIMPs and to learn more about their distribution in the galaxy. 
In addition, detection using different target nuclei would potentially al- 
low determination of both the spin-independent and spin-dependent cross 
sections. A measure of the direction of the recoiling nucleus would provide 
additional information on the distribution of WIMPs in the galaxy, allowing 
WIMP astronomy. 

Ultimately, the combination of WIMP direct and indirect detection 
with studies at colliders and of the cosmic microwave background could 
answer fundamental questions beyond whether WIMPs are the dark mat- 
ter and neutralinos exist. These combinations can determine whether the 
WIMPs are stable and whether there is non-baryonic dark matter other 
than WIMPs. As thermal relics, the WIMPs could provide a window to the 
early universe, or we may learn that the WIMPs must have been generated 
out of thermal equilibrium. WIMP astronomy could teach us about galaxy 
formation. Furthermore, the combination of information from these several 
methods may constrain the particle properties significantly better than col- 
liders alone. ^^'^^^ With the LHC taking data and more sensitive indirect 
and detect detection experiments in operation, the next five years may see 
an answer to the fundamental mystery of the nature of dark matter, as well 
as to other fundamental questions about the universe. 
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